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We previously developed two potent chemical classes that in- 
hibit the essential papain-like protease (PLpro) of severe acute 
respiratory syndrome coronavirus. In this study, we applied 
a novel approach to identify small fragments that act synergis- 
tically with these inhibitors. A fragment library was screened in 
combination with four previously developed lead inhibitors by 
fluorescence-based enzymatic assays. Several fragment com- 
pounds synergistically enhanced the inhibitory activity of the 
lead inhibitors by approximately an order of magnitude. Sur- 
face plasmon resonance measurements showed that three 


Introduction 


Human severe acute respiratory syndrome coronavirus (SARS- 
CoV) was initially reported in Guangdong province, China in 
2002, and quickly spread to almost 30 countries, resulting in 
about 8000 infections with approximately 800 deaths." It was 
contained through public health quarantine measures after 
several months. The isolation of closely related strains from the 
horseshoe bat has been reported,” which suggests a potential 
for reemergence of the SARS-CoV by retransmission to humans 
from its animal reservoir. There have also been three additional 
human coronaviruses—HCoV-HKU1,"! HCoV-NL63,"! and HCoV- 
EMC”!—reported since 2003. HCoV-HKU1 and HCoV-NL63 were 
shown to be less lethal than SARS-CoV.2* However, since 
April 2012, fifteen people have been infected with the most re- 
cently identified coronavirus, HCoV-EMC, nine of whom died 
from pneumonia and renal failure.” The whole genome se- 
quence of HCoV-EMC has been characterized, and it shows 
close similarity to two Asian bat coronaviruses: BtCoV-HKU4 
and BtCoV-HKU5.” This indicates a significant possibility for 
the emergence of new strains of SARS or new SARS-like 
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fragments bind specifically to the PLpro enzyme. Mode of in- 
hibition, computational solvent mapping, and molecular dock- 
ing studies suggest that these fragments bind adjacent to the 
binding site of the lead inhibitors and further stabilize the in- 
hibitor-bound state. We propose potential next-generation 
compounds based on a computational fragment-merging ap- 
proach. This approach provides an alternative strategy for lead 
optimization for cases in which direct co-crystallization is diffi- 
cult. 


human coronaviruses that could lead to even more deadly out- 
breaks. 

Human SARS-CoV has a 30-kb positive RNA genome and 
contains two proteases: 3-chymotrypsin-like protease (3CLpro) 
and papain-like protease (PLpro). PLpro and 3CLpro are re- 
sponsible for cleavage of the polyprotein into 16 nonstructural 
proteins. The first three positions between nonstructural pro- 
teins 1-4 are cleaved by PLpro, whereas the other eleven posi- 
tions are cleaved by 3CLpro. Both 3CLpro and PLpro have 
been shown to be essential for viral replication,”! making them 
attractive drug targets against SARS-CoV. 3CLpro has been in- 
tensively pursued as a drug target, and several inhibitors have 
been reported. They can be categorized as covalent peptido- 
mimetic, noncovalent peptidomimetic, and nonpeptidic."” 
Unlike 3CLpro, the development of PLpro inhibitors has been 
very limited until recently, despite its essential role in SARS 
viral replication. In addition to our own work,!'" two very differ- 
ent types of PLpro inhibitors have been reported in recent 
years. First, 6-mercaptopurine (6MP) and 6-thioguanine (6TG) 
were identified from compound screening, with IC,, values of 
5-20 um.|"7 Second, several natural products, diarylheptanoids 
isolated from Alnus japonica, showed inhibitory activity (IC., 
4.1 um) against PLpro."*! In addition to its primary function of 
viral peptide cleavage, PLpro has been recognized to be in- 
volved in deubiquitination, de-ISGylation, and viral evasion of 
the innate immune response.”*'“! SARS-PLpro is a cysteine pro- 
tease that contains a zinc binding motif, a catalytic triad, and 
a ubiquitin-like N-terminal domain (Figure 1). The sequence 
identity and similarity between SARS-PLpro and HCoV-EMC 
PLpro are 32 and 51%, respectively (Supporting Information 
(SI) figure S1a). Interestingly, the predicted structure of HCoV- 
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Figure 1. Overlay of the two known PLpro co-crystal structures. The inhibitor 2 bound (PDB ID: 3E9S) structure is 
shown in tan, and the inhibitor 3 bound (PDB ID: 3MJ5) structure is shown in cyan. The PLpro enzyme contains 
an intact zinc binding motif (left), a catalytic triad (middle), and a ubiquitin-like N-terminal domain. Inhibitor 2 and 
3 interactions with the PLpro enzyme are shown in the expanded square boxes in detail. 


EMC PLpro is very similar to that of SARS-PLpro (SI figure $1 b). 
The three catalytic residues (C112, H273, and D287) and four 
cysteines in the zinc binding motif are identical in these two 
PLpro enzymes. In addition, nine highly conserved residues in 
many other PLpro enzymes are also observed in HCoV-EMC 
PLpro (highlighted in yellow in SI figure S1a). This suggests 
that HCoV-EMC PLpro may possess characteristics similar to 
that of SARS-PLpro. 

By using high-throughput screening (HTS) and structure- 
based drug design, we previously developed a series of nonco- 
valent PLpro inhibitors that show promise as potential thera- 
peutic agents." Inhibitor 1 (GRL-0068S) has shown an inhibi- 
tory activity (IC,, value) of 0.29 um and antiviral activity (EC.9) 
of 5.2 um (see Table 1 below for inhibitor structures). Inhibitors 
2 (GRL-0617S), 3 (GRL-0667S), and 4 (GRL-0737S) have exhibit- 
ed inhibitory activities of 1.32, 0.64, and 1.36 um, with antiviral 
activities of 15, 9.1, and 9.1 um, respectively.'"'*'*! These inhibi- 
tors represent two active chemical scaffolds, with 1 and 2 rep- 
resenting the first (scaffold A), and 3 and 4 representing the 
second (scaffold B). Because these four lead inhibitors show 
moderate antiviral activity against SARS-CoV, their potencies 
must be further improved before development into clinical 
candidates. Analysis of the two available PLpro structures with 
inhibitors 2 and 3 and subsequent molecular modeling studies 
revealed that these leading inhibitors bind to a pocket formed 
mostly by the G267—G272 loop and the Leu163 side chain, dis- 
tinct from the catalytic site shown in green in Figure 1.'"'*! 
The inhibitors appear to block substrate access to the entrance 
of the catalytic site of the PLpro enzyme. Both scaffolds bind 
to a similar area, although scaffold A (inhibitors 1 and 2) 
blocks only a little over half of the area blocked by scaffold B 
(inhibitors 3 and 4). 

One way to improve the potency of these lead inhibitors 
would be to optimize the compounds by expanding into un- 
occupied space near their current binding site. This would 
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require the synthesis of larger 
compounds, a process that con- 
sumes_ significant amounts of 
time and resources, and may or 
may not result in better binding 
affinity and potency. An alterna- 
tive strategy is the application of 
fragment-based drug design 
(FBDD) methodology, which has 
become a_ well-established 
method for lead optimization." 
In many cases, FBDD has been 
used to discover small scaffolds 
for further optimization or to 
find two small fragments that 
can be linked together to pro- 
duce a more potent inhibitor.!' 
To be able to link two small frag- 
ments, it is typically necessary to 
solve crystal structures with each 
fragment. In many cases, obtain- 
ing this structural information 
can be very challenging and time-consuming due to the typi- 
cal weak binding affinity of fragments for their target enzymes. 
In this study, we used FBDD in combination with detailed en- 
zymatic analysis as a means to further optimize our leading in- 
hibitors without the use of crystallography. By testing a small 
library of fragment compounds in the presence of each of the 
lead inhibitors, we were able to identify fragments that synerg- 
istically improve the potency of the leading inhibitors. Addi- 
tional enzymatic analyses (including mechanism of inhibition 
and mutual exclusivity studies, described in detail below) and 
computational studies were used to predict the fragment bind- 
ing locations. Finally, a fragment-merging strategy was em- 
ployed to predict next-generation compounds that are able to 
recapitulate the binding position of the lead inhibitors as well 
as key binding features of the identified fragments. 


Ubiquitin-like domain 


Results and Discussion 
Fragment library screening 


The two currently developed SARS-CoV PLpro inhibitor scaf- 
folds do not bind at the catalytic site of the PLpro protease, 
but at a separate nearby binding site cradled between the 
L163 residue and binding loop 2 (G267—G272). This leaves the 
catalytic site unoccupied during inhibition. Although the cata- 
lytic pocket is small, there is a reasonable probability that frag- 
ment compounds could successfully bind to this location and 
enhance inhibitory activity of the lead inhibitors. The Zenobia 
library, consisting of 352 chemically and structurally diverse 
fragment-like compounds, was experimentally screened in the 
presence of each of the four previously discovered lead inhibi- 
tors. All fragment compounds were tested by continuous ki- 
netic enzymatic assays as described in the Experimental Sec- 
tion below. A replicate plot of inhibitor 2 with the 352 frag- 
ments showed highly reproducible results (SI figure S2a). The 
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replicate plots of the other lead inhibitors were similar to that 
of inhibitor 2. Of the 352 fragment-like compounds, 17, 14, 16, 
and 16 fragments showed > 30% increase in enzyme inhibi- 
tion when tested at a fragment concentration of 200 um in the 
presence of inhibitors 1, 2, 3, and 4, respectively. Percentage 
efficiency indices (PEI) for all 352 fragments were calculated as 
a hit selection criteria according to the method of Abad-Zapa- 
tero etal.; percent inhibition was divided by the molecular 
weight (M,) of each compound."'” Compounds with PEI values 
of 1.5 or greater were classified as hits, and the medium-range 
compounds (PEI 2.0-4.0) were selected for further testing. 
High PEl-range compounds (>4.0) were classified as outliers 
due to the strong likelihood of nonspecificity. The histogram of 
PEI values for all 352 fragments in combination with each of 
our four lead inhibitors is shown in SI figure S2b. A total of 13 
fragments were in the medium PEI range and were selected 
for further studies. 

The IC. values of the 13 selected fragments in the absence 
of any lead inhibitors were first determined. The IC,. values 
varied from 149 um (fragment F2) to 427 um (fragment F6), 
and are summarized in Table 1. Next, the IC., values of the four 


Table 1. 


1 (GRL-0068S) 


ant 


0.29 + 0.07 


Fragment IC, [UM] 1+ Frag. FD" 2-+Frag. 


F1 (ZT0009) 


188 +36 0.23 + 0.04 


2.94 + 0.09 


F2 (ZT0273) 


149+15 0.056 + 0.012 


0.66 + 0.01 


F3 (ZT0408) 


392 +28 0.049 + 0.021 


2.68 + 0.52 


F4 (ZT0426) 


212+27 0.032 + 0.013 


0.42 + 0.09 


F5 (ZT0470) 


388+10 0.050 + 0.011 


1.25 + 0.01 


F6 (ZT0537) 


427 +20 0.051 + 0.005 


0.53 + 0.06 


F7 (ZT0626) 


392-46 0.048 + 0.005 


0.34+0.13 


NH 
a. 0.036 + 0.013 


F8 (ZT0834) 


379 +38 0.43 + 0.08 


[a] Fold decrease in IC., value. 
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2 (GRL-0617S) 


es vouaen 


1.32 + 0.08 


lead inhibitors were determined in the presence of each of the 
13. fragment-like compounds. Of the 13 fragment-like com- 
pounds, eight fragments showed synergistic enhancement, 
with IC., value decreases of more than threefold for at least 
one of the lead inhibitors. The IC.. values of inhibitor 1 were 
enhanced 9.1- and 8.1-fold by fragments F4 and F8, respective- 
ly, relative to that of each lead inhibitor alone (Table 1). Frag- 
ment F7 showed the strongest synergy with inhibitor 2, en- 
hancing its IC,, value by 3.9-fold, and two other fragments (F4 
and F8) also equally enhanced inhibitor 2 by 3.1-fold. The IC. 
values of inhibitors 3 and 4 were enhanced 11.6- and 8.5-fold, 
respectively, by fragment F3, while the same compound, F3, 
did not improve the IC., value of inhibitor 2. Fragment F3 ap- 
peared to have better synergy with scaffold B inhibitors than 
with those representing scaffold A. Fragments F8 and F6 also 
showed significant synergistic activity with inhibitors 3 and 4, 
respectively. 

Six of the eight fragments exhibited reversible enzyme re- 
covery, whereas two—F3 and F5—showed respective recovery 
rates of only 25 and 31%. PLpro is a cysteine protease that 
contains a cysteine in the active site, and these partially irrever- 


IC,, value comparison of four lead inhibitors in combination with fragment compounds from the Zenobia library. 


Lead Inhibitor IC., [um] 


3 (GRL- oe 4 (GRL-0737S) 


0.64 + 0.06 
Fp! 


1.36 + 0.08 


FD" 4+ Frag. FDe 


0.20 + 0.04 0.79+0.19 


0.18 + 0.01 0.49 + 0.06 


0.055 + 0.027 0.16 + 0.03 


0.13 + 0.03 0.25 + 0.04 


0.17 + 0.04 0.44 + 0.08 


0.088 + 0.018 0.21 0.05 


0.12+0.01 0.26 + 0.03 


3.1 0.074 + 0.019 0.37 0.11 
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sible fragments may covalently a) 

, ; ; d 140 
interact with this residue. Cova- 

lently interacting compounds ue 
are likely to be nonspecifically 100 
reactive and usually have toxicity > 80 
o~ 
issues. Therefore, fragment F3 < 60 


may not be a good candidate, 40 
although it was the best enhanc- 
er for scaffold B inhibitors (inhib- 
itors 3 and 4). Similarly, these re- 
sults indicate that fragment F5 0 
may also exhibit problems with 
irreversibility. Two — ligand-effi- 
ciency-related indices, binding 
efficiency index (BEl) and surface 
binding efficiency index (SE)), 
were calculated as described.""®! 
BEI relates potency to the molec- 
ular weight of each compound, 
while SEI is an index that moni- 
tors the potency gains as polar 
surface area (PSA) increases. BEI 


and SEI for these eight frag- ) 200 400 
[Fragment F3] / um 


ments and four lead inhibitors 
were calculated by simple equa- 
tions: pIC., divided by M, for BEI 
and plC., divided by PSA for 
SEI."'7 The PSA value is an impor- 
tant indicator related to oral bio- 
availability and intestinal perme- 
ability." The BEI values of fragments F3, F4, and F7 were > 25, 
suggesting that they are outliers due to the possibility of non- 
specific interaction (SI figure S2c). Fragment F7 contains an al- 
dehyde, which is likely to exhibit nonspecific affinity for the 
cysteine residue in the catalytic site of the PLpro enzyme. 
Thus, considering both ligand efficiency indices and enzyme 
activity recovery rate, fragments F8, F8, F6/F8, and F6 are the 
best enhancers for inhibitors 1, 2, 3, and 4, respectively. 


tails). 


Binding analysis by SPR 


Although we observed IC., value enhancement by the frag- 
ments through dose-response curves using enzymatic assays, 
it is possible that the results could be due to nonspecific inter- 
actions between the PLpro enzyme and the _ fragments. 
Enzyme assays alone can lead to false positives in some cases. 
To investigate this further, orthogonal binding analyses by sur- 
face plasmon resonance (SPR) were performed, and the disso- 
ciation equilibrium constants (Kp), which reflect binding affinity, 
were determined for each of the four lead inhibitors 1-4 and 
the seven newly identified fragments F2-F8. Fragment F1 was 
excluded from these studies due to its low IC., enhancement. 
The binding affinities of inhibitors 1-4 varied between 0.97 
and 1.67 um as determined by SPR, which were in a range sim- 
ilar to their IC., values (Table 2 and Figure 2a). Promisingly, the 
binding affinities of fragments F2, F6, and F8 were respectively 
202, 461, and 437 um. The K, fitting curve of F8 is shown in 
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[Inhibitor 2] / uM 
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[Fragment F8] / uM 
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[Fragment F5] / um 


Figure 2. Dissociation equilibrium constant (Kp) determination by SPR binding analysis: fitting curves of a) inhibi- 
tor 2 and b) fragment F8. The response unit difference (ARU) plots of c) fragments F3 and d) F5 with respect to 
concentration of each compound are shown, and suggest random and nonspecific binding, respectively. Data in 
panels a) and b) were fit to a single rectangular hyperbolic curve to calculate Kp (see Experimental Section for de- 


Table 2. Dissociation equilibrium constants (Kp determined by SPR) and 
mechanism of inhibition. 


Compound Kp [um] Inhib. Type K; [um] 


1 (GRL-0068S 
2 (GRL-0617S 
3 (GRL-0667S 
4 (GRL-0737S 
F2 (ZT0273) 
F6 (ZT0537) 
F8 (ZT0834) 


1.69 + 0.27 
1.04 + 0.16 
0.97 + 0.03 
1.06 + 0.09 

202 + 89 

461 + 30 

437 = 36 


noncompetitive (a = 2.6 
noncompetitive (a = 2.0 
noncompetitive ( 
noncompetitive ( 
noncompetitive ( 
noncompetitive ( 
noncompetitive ( 


a= 
a= 
a= 
a= 
a= 


Figure 2b as an example. However, three fragments, F3, F4, 
and F7, were suspected of forming nonspecific interactions 
due to their unusually high BEI values as noted above. SPR 
analyses confirmed that these three fragments, as well as frag- 
ment F5, all bind nonspecifically. Fragments F3, F4, and F7 did 
not show a typical dose-response curve as their concentra- 
tions increased; instead, they showed a random binding pat- 
tern as shown Figure 2c for F3. Fragment F5 did show a typical 
dose-response curve, but SPR response units increased to 
nearly 4000 RU when its concentration increased to 800 um 
(Figure 2d). This is likely due to multisite binding caused by 
aggregation rather than one-to-one interaction with the PLpro 
enzyme, suggesting that the fragment may bind to many dif- 
ferent locations on the PLpro enzyme. 
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Mechanism of inhibition studies 


Mechanism of inhibition studies have not been thoroughly de- 
scribed for the lead inhibitors. Prior studies suggested that in- 
hibitors 2 and 3 behave as competitive inhibitors."*'! Two 
crystal structures with these compounds verify that they bind 
at the entrance of the catalytic site, blocking substrate access 
to the site. Although these lead inhibitors do not bind directly 
to the catalytic site, they nonetheless display competitive in- 
hibition patterns. To thoroughly investigate the inhibition 
mechanism of both lead inhibitors and fragment compounds, 
kinetic studies for each compound were performed with the 
enzyme-inhibitor complexes and varying substrate concentra- 
tions. The data were fitted to four equations [Eqs. (2)—(5) in Ex- 
perimental Section] and three plots (Michaelis-Menten, Line- 
weaver-Burk, and Dixon), each using SigmaPlot Enzyme kinet- 
ics Module 1.3. Akaike information criterion-corrected (AICc) 
values were used to determine the best fit equation?” for 
which the lowest AlCc value corresponds to the best fit equa- 
tion. To determine which equation fits best, the AlCc value dif- 
ference must be at least two units from the next lowest.” All 
four previously developed lead inhibitors showed noncompeti- 
tive inhibition with unequal binding affinity to free enzyme (E) 
and substrate-bound enzyme (ES), where the alpha (qa) value is 
not equal to 1 (Table 2). All four lead inhibitors showed 
a values greater than 1, which indicates the binding affinity to 
free enzyme is tighter than that of ES. Dixon plots for noncom- 
petitive inhibition of inhibitor 1 (a=2.6) with respect to the 
substrate are shown in Figure 3a, and the resulting Kk; value 
was 0.42 um. All three fragments showed noncompetitive in- 
hibition with different a values as well. Fragments F6 and F8 
had a values greater than 1 (a> 1), indicating that these frag- 
ments also behave similarly to the lead inhibitors, although 
a values of these fragments were smaller than those of the 
lead inhibitors. This may explain why F8 displayed the stron- 
gest effect with the three lead inhibitors (1, 2, and 3), whereas 
fragment F6 exhibited the strongest synergistic effect with 
scaffold B inhibitors (3 and 4). Dixon plots for noncompetitive 
inhibition of fragment F8 (a= 1.6) are shown in Figure 3b, and 
the K, value was determined to be 242 um. Because all four 
lead inhibitors and two fragments showed noncompetitive in- 
hibition with a>1, it is likely that these compounds bind to 
another site that is not a catalytic site, hindering substrate 
binding to the enzyme. Therefore, these studies suggest that it 
is unlikely that these newly identified fragment-like com- 
pounds bind to the catalytic site of the PLpro enzyme. 


Computational solvent mapping 


From our mechanism of inhibition and binding synergy/ 
mutual exclusivity analyses, we determined that the newly 
identified fragment-like compounds bind to a site separate 
from that of the lead inhibitors. Unfortunately this information 
is insufficient to identify where the fragments bind on the 
PLpro enzyme. Therefore, we investigated all possible binding 
site candidates. The unexplored catalytic pocket was designat- 
ed as candidate fragment binding site 1, although our mecha- 
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Figure 3. Inhibition mode of lead inhibitors and fragments: Dixon plots for 
noncompetitive inhibition of a) inhibitor 1 (GRL-0068) and b) fragment F8 
with respect to the substrate Z-Arg-Leu-Arg-Gly-Gly-AMC at indicated con- 
centrations. K; values determined for inhibitor 1 and fragment F8 were 0.42 
and 241 um, respectively. Equations (2)-(5) (see the Experimental Section 
section) in the SigmaPlot Enzyme Kinetics Module 1.3 were used to fit the 
experimental data. The noncompetitive inhibition model was the best fit for 
both inhibitor 1 and fragment F8, producing respective a values of 2.6 and 
1.6. 


nism of inhibition studies indicated that this location was not 
likely to be the fragment binding site. To identify additional 
candidate fragment binding sites, we performed a series of 
computational solvent mapping experiments, or “hot spot” 
analyses, using the FTMAP server (see Experimental Section).?" 
FTMAP is a multistage protein mapping algorithm that is 
based on a fast Fourier transform (FFT) correlation. This ap- 
proach can efficiently search for potential binding sites on the 
entire surface of the protein. Two model systems, crystal struc- 
tures of inhibitor 2 bound to PLpro (PDB ID: 3E9S) and inhibi- 
tor 3 bound to PLpro (PDB ID: 3MJ5), were used for the analy- 
ses. These two provided representative crystal structures for 
each compound scaffold. Using 50ns molecular dynamics 
(MD) simulations for each system (3E9S, 3MJ5), we determined 
that the major fluctuations of both structures come from the 
zinc binding motif and a B sheet forming the catalytic triad. 
Representative snapshot structures were extracted from the 


ChemMedChem 2013, 8, 1361-1372 1365 


MED 


simulations as discussed below (see Experimental Section) and 
submitted to the FTMAP server, which performed a fragment- 
based binding site analysis. Two strong candidate fragment 
binding sites were identified by analyzing the original struc- 
tures and consensus clusters of the probe molecules used (see 
Experimental Section below). The two hot spots identified by 
this analysis included an extension of the original binding site 
of the lead inhibitors and one cavity in the palm region 
(shown in Figure 4a). Hot spot 1 was larger than the lead in- 
hibitor binding site and included an additional cavity unoccu- 
pied by the lead inhibitors. This binding site has been previ- 
ously discussed as a possible substrate recognition site.' This 
extended cavity is composed of five residues: R167, E168, 
M209, D303, and T302. The second position identified by 
FTMAP, hot spot 2, was located >10A away from both the 
lead inhibitor binding site and the catalytic site of PLpro. It is 
encompassed by the zinc binding motif and the palm region. 
Lastly, the small volume containing the catalytic triad (site 1) 


a) 


Hot spot 1 


; Hot spot 2 


Binding Site 2 
‘ (Hot spot 1) 


b) C) 


Binding Site 1 Binding Site 3 
(Catalytic site) (Coarse docking) 


Figure 4. Fragment binding site analysis results from FTMAP and identifica- 
tion of potential binding sites: a) Shown are two potential binding sites 
identified by FTMAP using the co-crystal structure of PLpro with inhibitor 3. 
The first candidate site (upper circle) is an extension of the lead inhibitor’s 
binding site, formed by E186 (grey), R187 (orange), M209 (blue), T302 
(green), and D303 (pink). The second candidate site (lower circle) is encom- 
passed by the zinc binding motif and the palm region. b) Unoccupied cata- 
lytic pocket surrounded by three residues (C112, H273, and D287). This small 
pocket was initially considered a potential fragment binding site 1, but sub- 
sequently eliminated. c) Two small unoccupied pockets near the PLpro lead 
inhibitor binding site. Binding site 2 was identified by FTMAP hot spot analy- 
sis, and the third putative binding site (binding site 3) was identified by 
docking studies in the region near the lead inhibitor binding site. The image 
in panel c) was prepared by 90° counterclockwise rotation of the image in 
panel b), which gives a better view of binding site 2. 
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was not identified by FTMAP, likely due to the smaller size of 
this pocket. 


Fragment binding site analysis 


Structural analysis and computational studies identified four 
potential fragment binding sites, three of which are located 
within 10 A of the binding sites of the lead inhibitors. Candi- 
date binding site 1, shown in green, is the unoccupied catalytic 
pocket (Figure 4b), identified by analysis of the co-crystal 
structures available. Candidate binding site 2 was identified by 
FTMAP using computational solvent mapping, as discussed 
above and in the Experimental Section. The surrounding resi- 
dues of candidate binding site 2, located to the right of the 
lead inhibitor binding site in Figure 4c, are shown in blue. Ad- 
ditional coarse docking studies with hit fragment compounds 
(see Experimental Section) identified another shallow binding 
site (designated candidate binding site 3, shown in orange in 
Figure 4c). The mechanism of inhibition analyses and mutual 
exclusivity studies discussed above indicated that the frag- 
ments bind in combination with the lead inhibitors in a distinct 
and separate binding site of the enzyme. This important obser- 
vation lays the foundation for future fragment-inhibitor merg- 
ing or linking studies proposed below. 

To further investigate the candidate fragment binding sites 
and to predict the most likely binding conformation of the 
fragments within each site, we performed a series of computa- 
tional binding free energy calculations. These studies, de- 
scribed in detail in the Experimental Section, used a combina- 
tion of molecular docking, molecular simulations, and molecu- 
lar mechanics/Poisson—Boltzmann surface area (MM/PBSA) cal- 
culations to predict the binding free energy (AG,,,4) of the 
lead inhibitors when the synergistic fragments were placed in 
the adjacent candidate binding sites. These studies were per- 
formed using three potential binding sites: site 1 (the catalytic 
site), site 2 (the extended lead inhibitor binding site), and site 3 
identified by coarse docking. Hot spot 2 was excluded from 
these studies due to the large distance (> 10 A) from the inhib- 
itor binding site. Hot spot 2 is very shallow, highly solvent- 
exposed, and located between the zinc binding motif and the 
palm region (see Figure 4a). More importantly, it appears likely 
that even if binding of small fragments to this position could 
cause some degree of conformational changes to these two 
regions, any structural change would not be likely to affect the 
lead inhibitor binding site to enhance their IC,, values to the 
extent observed. For these reasons, as well as the impracticali- 
ty of fragment linking or merging strategies at such a distance, 
hot spot 2 was excluded from further analyses. 

The calculated free energies of binding of the lead inhibi- 
tors, with a fragment placed in each successive candidate site, 
were then compared with the experimentally determined bind- 
ing affinity (IC..) of the lead inhibitors with the synergistic frag- 
ments to predict the most likely fragment binding site. For the 
purposes of these studies, the binding conformations of inhibi- 
tors 2 and 3 from the co-crystal structures, representing scaf- 
folds A and B, were used. Fragment F8 was selected for use in 
these studies due to its strong synergistic effect on both inhib- 
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itors 2 and 3. The fragment was docked into each of the three 
candidate binding sites discussed above, and three docking 
poses for each fragment were selected for further simulation 
(see Experimental Section). Post-analysis of the simulation re- 
sults by MM/PBSA allowed calculation of the binding free ener- 
gies of the lead inhibitors for each simulation, as shown in 
Figure 5. 


AGping / kcal mol-' 


C1 C2 C3 C1 C2 C3 C1 C2 C3 C1 C2 C3 
Binding site 1 Binding site 2 Binding site 1 Binding site 2 


Inhibitor 2 + F8 Inhibitor 3 + F8 


Figure 5. Calculated binding free energies of inhibitors 2 and 3 when frag- 
ment F8 is docked into various candidate sites. Three different docked con- 
formations, C1 (circles), C2 (triangles), and C3 (squares) at each binding site 
were selected for calculations. More negative numbers represent better 
binding conformation, with those marked by black arrows being the best. 


We initially noted that fragment F8 was not stable in bind- 
ing site 3, and repeatedly drifted from the site by the 20 ns 
mark, regardless of its starting position. This suggests that the 
shallow and solvent-exposed binding site 3 is also not an ap- 
propriate position for fragment binding. Hence, the site was 
also removed from consideration, and further analysis of this 
site was not performed. Additionally, an analysis of the MM/ 
PBSA results (Figure 5, SI tables $1 and S2) revealed that the 
catalytic site, site 1, was also not a favorable position for frag- 
ment binding. The binding free energy values for the inhibitor 
with fragment F8 at this position were higher than those cal- 
culated for the inhibitor alone as well as the experimentally de- 
rived results. These results, taken together with the experimen- 
tal results discussed above, indicate that the catalytic site is 
most likely not the binding site of the identified synergistic 
fragments. Favorable results that closely aligned with experi- 
mental values were obtained, however, when the fragments 
were placed at binding site 2 (shown in Figure 6), consistent 
with synergistic activity. We noted a preference for a positive 
charge, as observed in the F8 piperazine amine group, to be 
placed into an anionic pocket of site 2, defined by D303, T302, 
Y274, and D165, which, in combination with the established 
binding features of the lead inhibitors, resulted in a clear gain 
in their calculated binding affinity. An additional edge-to-face 
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Figure 6. The most probable binding conformation of fragment F8 (orange) 
in combination with lead inhibitors a) 2 (magenta) and b) 3 (magenta). Resi- 
dues of the PLpro enzyme that interact with fragment F8 are shown in the 
boxes at right, with hydrogen bonds shown in green. 


m stacking interaction between the naphthalene ring system 
and the aromatic ring system of fragment F8 (visible in 
Figure 6) may further stabilize the lead inhibitors in their bind- 
ing conformations. 

To prioritize compound candidates for future synthesis and 
testing, we undertook a series of fragment-merging and dock- 
ing analyses, using the information determined from the bind- 
ing site prediction studies discussed above. Our fragment- 
merging strategy involved the elimination of a potentially 
problematic nitroaromatic ring of fragment F8 and the direct 
incorporation of the piperazine ring into the lead inhibitor 
scaffold, as shown in Figure 7. A series of proposed com- 
pounds were designed using a variety of linking groups and 
positions (discussed in the Experimental Section below), and 
docked into the extended binding site defined by the lead in- 
hibitor and binding site 2 (Table 3). Proposed compounds were 
selected by score and their ability to successfully recapitulate 
the key binding features of the inhibitor and synergistic frag- 
ment, including the lipophilic interactions of the naphthalene 
ring system and the positioning of the positively charged 
amine of the piperazine into the key anionic pocket discussed 
above. 


Conclusions 


In attempts to improve the potency of currently developed 
leading PLpro inhibitors, a fragment-like library was enzymati- 
cally tested in combination with four lead inhibitors. We dis- 
covered a total of eight fragment-like compounds that synerg- 
istically enhanced the inhibitory activity of lead inhibitors. An 
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Figure 7. One of the compounds proposed for synthesis and testing by the 
fragment-merging strategy: The proposed compound M1 (green carbon 
atoms), after molecular docking, is able to recapitulate the crystallographic 
position of the inhibitor and the predicted position of the synergistic frag- 
ment (orange carbons). The potentially problematic nitroaromatic ring of 
the fragment is removed in favor of the merged compound, which directly 
incorporates the piperazine ring into the lead inhibitor scaffold. Key binding 
features of the inhibitor and synergistic fragment are conserved, including 
the lipophilic interactions of the naphthalene ring system and the position- 
ing of the positively charged amine of the piperazine into a key anionic 
pocket. 


Table 3. Compounds proposed for synthesis and testing by fragment 
merging and docking calculations. 


GRL-0617S + ZT0834 


Compd Structure Docking Score 


[a] Structure portions not highlighted are from the lead inhibitors, and 
moieties highlighted in grey circles came from fragment F8. Linkers be- 
tween the lead inhibitor and a fragment part are highlighted with a small 
grey oval. 
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orthogonal binding analysis by SPR confirmed the binding of 
three of these fragments and also determined their direct 
binding affinities (Kp). Our characterization assays of the bind- 
ing mechanism of the two newly identified fragment com- 
pounds showed noncompetitive inhibition. This indicated that 
the synergistic fragment compounds did not bind to the cata- 
lytic site of the PLpro enzyme. Computational analysis predict- 
ed several potential fragment binding sites, and more expan- 
sive studies revealed a preference for fragment binding to 
site 2, an extension of the lead inhibitor binding site, when the 
inhibitor binding site was occupied by representative struc- 
tures. Enzymatic characterization and computational studies 
suggest that these fragments, by binding adjacent to the bind- 
ing site of the lead inhibitors, are able to further stabilize the 
inhibitor-bound state. Using the structural information § re- 
vealed by these studies, we were able to apply a fragment- 
merging strategy to suggest several next-generation com- 
pounds with high predicted binding affinity. The novel ap- 
proach described herein may be beneficial for other cases in 
which direct crystallization studies prove difficult. 


Experimental Section 


Reagents and chemicals: The Zenobia library, consisting of 352 
structurally diverse small fragment-like compounds, was purchased 
from Zenobia Therapeutics. The molecular weights of these 352 
compounds vary between 94.12 and 285.9 Da, with a mean molec- 
ular weight of 154Da. The mean number of hydrogen bond 
donors and acceptors are 1.4 and 2.6, respectively. The mean 
clogP and topological polar surface area (tPSA) values are 1.6 and 
52 A’, respectively. The solubility of each compound was verified 
by the vendor to be at least 200 mm in 100% DMSO. The 352 com- 
pounds in the Zenobia library are functionally and chemically di- 
verse and synthetically accessible for further optimization. They are 
categorized and distributed as shape-diverse mixtures of eight 
compounds in 96-well plate columns. The four lead inhibitors 1-4 
were provided by Dr. A. K. Ghosh. The PLpro substrate, Z-Arg-Leu- 
Arg-Gly-Gly-AMC, was purchased from Bachem Bioscience. 


Tested compounds: Previously developed inhibitors 1-4 were syn- 
thesized by Dr. A. K. Ghosh, and the purity of each compound was 
confirmed by NMR as described.''”'*! Fragments F1-F8 were re- 
purchased from Sigma. Compound purity was determined by NMR 
to be >95%. NMR spectra for each fragment compound listed in 
Table 1 are included in the Supporting Information. 


Cloning, expression, and purification of SARS-CoV PLpro: The 
PLpro gene (SARS-CoV polyprotein residues 1541-1855) with 
human rhinovirus 3C protease cleavage site and His tag 
(LEVLFQGPHHHHHH) at the C terminus was initially prepared by 
codon-optimized gene synthesis (BioBasic Inc.) and cloned into 
a pET15b vector between Necol and Xhol sites. A stop codon was 
added immediately after the last PLpro gene to generate native 
PLpro without any tag. Rosetta2(DE3) cells (Novagen) containing 
the recombinant plasmid were grown to an OD, of 0.6 at 37°C 
by shaking at 200 rpm in LB medium (2 L volume) containing am- 
picillin (100 uwgmL~') and chloramphenicol (34 ugmL'). Native 
PLpro was expressed by adding 0.5 mm IPTG with incubation con- 
tinued at 25°C for 4h. The cells were then harvested by centrifu- 
gation at 8000 rpm (Beckman JA-10) for 15 min and resuspended 
in lysis buffer (1 mgmL7! lysozyme and protease inhibitor cocktail 
in buffer A: 20mm Tris‘HCl (pH 7.5), 10mm _ £-mercaptoethanol 
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(BME)) and were sonicated on ice. The lysate was centrifuged at 
19000 rom (Beckman JA-20) for 30 min, and the supernatant was 
collected. The supernatant was then subjected to a 40% ammoni- 
um sulfate fractionation. The suspension was centrifuged at 
19000 rom (Beckman JA-20) for 20 min, and the resulting pellet 
was resuspended in buffer B (1.5m ammonium sulfate, 20 mm Tris, 
pH 7.5, 10 mm BME) and centrifuged at 19000 rpm (Beckman JA- 
20) for 20 min to remove precipitate. The supernatant was loaded 
onto a 5 mL phenyl-Sepharose 6 Fast-Flow HS column (Amersham 
BioSciences, Piscataway, NJ, USA) equilibrated with buffer B. The 
protein was eluted with a stepwise gradient of 20, 40, 60, and 
100% buffer A. The native PLpro eluted at 100% buffer A. Fractions 
containing PLpro were pooled and dialyzed with buffer A over- 
night at 4°C and loaded onto a HiTrap QXL (Amersham BioSci- 
ences). The protein was eluted with a stepwise gradient with buf- 
ferC (20mm Tris, pH 7.5, 0.5m NaCl, 10mm BME), and PLpro 
eluted out at ~20% buffer C. Pure SARS-CoV PLpro fractions were 
combined, buffer-exchanged into buffer A containing 20% glycer- 
ol, and stored at —80°C. Purity was >90%, as determined by SDS- 
PAGE analysis. 


Compound screening: The Zenobia library was purchased as 
200 mm DMSO stocks. All 352 compounds were experimentally 
tested against PLpro in combination with the four previously de- 
veloped PLpro inhibitors. PLpro activity was measured by a continu- 
ous kinetic assay with ubiquitin-derived substrate Z-Arg-Leu-Arg- 
Gly-Gly-AMC (Bachem Bioscience) that generates a fluorescence 
signal (excitation A: 360 nm; emission 2: 450 nm) when aminome- 
thylcoumarin (AMC) is cleaved. For primary screening, Zenobia 
compounds were first diluted to 10 mm DMSO stock in a 384-well 
plate (Corning Inc.) from 200mm original compound solutions. 
Assay buffer (7.5 uL) containing 50 mm HEPES, pH 7.5, 0.01 % Triton 
X-100, 0.1 mgmL~' BSA, and 2mm GSH was distributed to black 
low-volume 384-well plates (Corning Inc.). Zenobia compounds 
(0.6 wL, 10 mm) were added to the assay buffer. PLoro enzyme so- 
lution (10 nm final concentration) was prepared in the presence of 
each of the four lead inhibitors at a concentration equal to its IC., 
value. PLpro in the absence of any inhibitors was also prepared as 
controls. Enzyme-lead inhibitor mixture solution (10 uL) was then 
added to the plates containing the Zenobia compounds, and the 
plates were incubated for an additional 5-10 min. The reaction 
was initiated by adding substrate (7.5 uL, 75 um final concentra- 
tion), and fluorescence intensity was continuously monitored with 
a POLARstar OPTIMA microplate reader (BMG LABTECH) for 10 min. 
All compounds were tested in duplicate, and each plate contained 
a total of 16 positive and 16 negative controls. 


IC, value determinations: The IC,, values of each fragment com- 
pound that showed synergistic inhibitory activity were determined. 
The hit compounds were prepared at various concentrations (0- 
1000 um). The IC,, values of the four lead inhibitors were measured 
in the same concentrations of enzyme and substrate as the initial 
screening with a series of the lead compound concentrations (0- 
20 um) in the absence and presence of 200 um hit compounds. 
The reaction was initiated by adding fluorogenic substrate, and its 
activity was continuously monitored. The IC.. values were calculat- 
ed by fitting to the three-parameter Hill equation [Eq. (1)] with Ori- 
ginPro 8.5 (OriginLab, Inc.): 


x” 
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for which y is percent inhibition, x is inhibitor concentration, n is 
the slope of the concentration-response curve (Hill slope), and 
Vinax IS Maximal inhibition from three independent assays. 


Reversibility of fragment inhibition: The reversibility of the hit 
fragment compounds was determined by the dilution method.” 
PLpro enzyme was prepared as a 400 nm solution (20x assay con- 
centration) and was incubated with screened compounds at 20x 
each IC., concentration for 30 min at room temperature in assay 
buffer containing 50 mm HEPES pH 7.5, 2mm GSH, 0.01% Triton, 
0.1 mgmL7'! BSA, and 1% DMSO in a final volume of 200 uL. PLpro 
with the same amount of DMSO in place of fragment compounds 
was also prepared in the same way as a control. The activity of 
PLpro was measured in the same way as IC., measurements. The 
enzyme-inhibitor solution was then diluted 20- and 40-fold and in- 
cubated for 30 min before measurement of percent recovery of 
enzyme activity. All reversibility assays were done in triplicate. 


Determination of dissociation equilibrium constant (K,) by SPR: 
Untagged native PLpro enzyme was prepared in PBS (10 mm phos- 
phate pH 7.4, 2.7mm KCI, 137mm NaCl) and immobilized on 
a CM5 sensor chip using standard amine coupling at 20°C with 
running buffer HBS-P (10 mm HEPES, 150 mm NaCl, 0.05% surfac- 
tant P-20, pH 7.4) from GE healthcare using a Biacore T100 instru- 
ment. Flow channels 1 and 3 were activated by a mixture of 1- 
ethyl-3-(3-dimethylaminopropyl)carbodiimide hydrochloride (EDC)/ 
N-hydroxy succinimide (NHS), and the activated surface was 
blocked by ethanolamine (pH 8.5) as controls. PLpro enzyme was 
diluted in 10 mm sodium acetate (pH 5.0) and immobilized to flow 
channels 2 and 4 after sensor surface activation with EDC/NHS 
with a 7 min injection followed by ethanolamine blocking on unoc- 
cupied surface area. PLpro immobilization levels of flow channels 2 
and 4 were ~12500 and ~8500 RU, respectively. The four previous- 
ly discovered inhibitors and seven identified Zenobia fragment 
compounds were prepared as 10 mm (inhibitors 1-4) or 200 mm 
(fragments) DMSO stocks. Compound solutions with a series of in- 
creasing concentrations (0-20 um for inhibitors, 0-800 um for frag- 
ments at 1.5-fold dilution) were applied to all four channels at 
a flow rate of 10 uLmin™'. Sensorgrams were analyzed using Biae- 
valuation software 2.0.3, and response unit difference (ARU) values 
at each concentration were measured during the binding equilibra- 
tion phase. SigmaPlot 11.0 was used to fit the data to a single rec- 
tangular hyperbolic curve to determine Kp values. The hyperbola 
function, Y=Yn’X/(Kp-+X), Was used to plot response units and cor- 
responding concentration, where y is the response, y,,,, is the max- 
imum response, and x is the compound concentration. 


Mechanism of inhibition: PLpro activity was monitored in the 
same way as the primary screening with varying concentrations of 
both inhibitor compounds and substrate (0-300 um). The concen- 
tration of compounds was varied from 0 to at least 10x the IC. 
value of each fragment. The data were fit to the following equa- 
tions using SigmaPlot Enzyme Kinetics Module 1.3 to determine 
the best fit inhibition mechanism and kinetic parameters for each 
compound: 


Competitive : | vy = ——_“"-_—__ 

14+2)+(1+_ 2) 
[S] [Ki] 
Noncompetitive Vax 
itive : vy = ———__"“"*___ 
3 
(+m) +(1+8) (3) 
e,e Viegas 
Uncompetitive : vy = ———_- (4) 
1+ K] - BE 
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Mixed-type: v= 


(i) (14m) + (+a) " 


for which v is the reaction rate, V,,,, is the maximum rate of the re- 
action, K,, is the Michaelis-Menten constant for substrate, [S] is the 
substrate concentration, [I] is the inhibitor concentration, K, is the 
dissociation constant of the inhibitor | from the free enzyme, and 
ak, is the dissociation constant for the inhibitor | from the ES com- 
plex. 


Computational studies 


Fragment binding site analysis: All MD simulations were per- 
formed using the AMBER suite of programs (AMBER 11).”7! The two 
known co-crystal structures (inhibitor-bound, PDB IDs: 3E9S and 
3MJ5) were obtained from the RCSB Protein Data Bank. Multimer 
protein complexes were reduced to the functional monomer form 
for simulations. The bound inhibitor and zinc ions were retained in 
the systems; all other ions and crystal waters were removed. The 
ligand parameters were generated using the Antechamber soft- 
ware in AmberTools v1.5 and the General Amber Force Field 
(GAFF) with am1bcc charges.” Leap was used to neutralize the 
systems by adding sodium ions, and each system was solvated 
using 8 A TIP3P water in a rectangular box. The systems were pre- 
pared using the ff99SB Force Field.” Minimization using SANDER 
was carried out in two steps: first, only the water molecules were 
minimized for 5000 steps while keeping the protein fixed; second, 
the entire system was minimized for an additional 10000 steps. Pe- 
triodic boundary conditions were applied with 8 A of nonbonded 
cutoff. The systems were heated from 0 to 300 K using the NVT en- 
semble in 50 ps with positional restraint on the protein. Following 
this, the systems were equilibrated for 100 ps using the NPT en- 
semble without restraint. The production phase of the simulation 
was performed using the NPT ensemble without positional re- 
straint for 50 ns. Post-analysis was performed using the _ ptraj 
module from AmberTools v1.5 to compute the average structures 
at every 1 ns and write out the PDB file for the averaged coordi- 
nates. A total of 50 snapshots were extracted using this method. 


The PDB files generated above were uploaded to the FTMAP 
server (http://ftmap.bu.edu/). The FTMAP server removes any 
waters, ions, and bound ligands and only retains the protein struc- 
ture. Sixteen different probes were used: ethanol, isopropanol, iso- 
butanol, acetone, acetaldehyde, dimethyl ether, cyclohexane, 
ethane, acetonitrile, urea, methylamine, phenol, benzaldehyde, 
benzene, acetamide, and N,N-dimethylformamide. The output re- 
sults generated by FTMAP were in the form of PDB files with con- 
sensus clusters of the different probes mapped to the surface of 
the PLpro structure. Cluster sites identified by FTMAP were re- 
quired to be conserved among all molecular probes for selection 
and further study, described below. Two such sites were identified 
(see Discussion above). One additional candidate binding site, 
site 3, was identified using coarse docking studies, described 
below. 


Binding affinity predictions: To further investigate potential bind- 
ing sites, fragments were docked to the enzyme using Glide-XP 
(extra precision) and the Induced Fit Docking (IFD) protocol of 
Schrédinger, using default settings." The structures of PLpro 
bound with GRL-0617S (PDB ID: 3E9S) and GRL-0667S (PDB ID: 
3MJ5) were obtained from the RCSB Protein Data Bank. Prior to 
molecular docking, structures of protein targets were prepared 
using the Protein Preparation Wizard in the Schrédinger software 


© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim 


package (Suite 2012). All water molecules in the crystal structures 
were removed, and crystal multimers were reduced to the mono- 
mer. Residues located within 5 A of the bound inhibitors were re- 
fined by the Prime algorithm. Two distinct docking runs were per- 
formed using this strategy. The first docking run was blind docking 
(coarse docking), with the docking area defined to be within 10 A 
of the lead inhibitors. Candidate site 3 was identified during these 
studies as a high-scoring binding site for multiple fragments. The 
subsequent docking runs were defined by the location of the can- 
didate sites identified by the structural and computational analy- 
ses, described above. Binding site 1 (the catalytic site) was defined 
by selecting residues H273, C112, and D287; binding site 2 was de- 
fined by selecting residues R167, E168, M209, D303, and T302; and 
binding site 3 was defined by selecting residues N157, K158, D165, 
and E168. Binding hot spot 2 was not studied for the reasons de- 
scribed above. In each case, the docking area was defined to be 
a 10A box centered on the selected residues. Fragments of the 
Zenobia library were prepared in multiple tautomer and protona- 
tion states at pH 7.4 by using LigPrep (Schrédinger).?7 All eight 
active fragments were docked in the coarse docking studies; only 
fragment F8 was docked in the subsequent candidate site docking 
studies. All docking poses were ranked by Glide XP score.”®! In the 
second docking run, three conformations (poses) each of F8 at 
each binding site (sites 1, 2, and 4) produced by IFD were selected 
as the starting conformations for the MD simulations and subse- 
quent MM/PBSA calculation, described below. The binding poses 
selected included the top score pose and two other poses which 
possessed the greatest difference (by calculated RMSD). 


MD simulations and subsequent MM/PBSA calculations used to 
predict binding affinities were performed using the AMBER suite of 
programs (AMBER 11).”*! Preparation of the MD simulations, includ- 
ing parameterization of the fragment (F8) and inhibitors, and prep- 
aration and equilibration of the system, were performed exactly as 
described above for binding site analyses. The production phase of 
the simulation was performed using the NPT ensemble without po- 
sitional restraints for 20 ns. The python script, MMPBSA.py, includ- 
ed in AMBER 11, was used to perform the MM/PBSA calculations.”” 
The first 100 frames from the last 2 ns (equilibrated) of the simula- 
tion trajectory file were used for the MM/PBSA calculations. The 
ionic strength was set to 0.1m, and the internal dielectric constant 
was set as 1.0 (default). All other options were set to default set- 
tings. The values for the free energy of binding of each inhibitor, 
combined with the three different fragment poses, were calculated 
according to Equation (6): 


AG ping = Ganipies— Greceptat= Ciaand (6) 


in which Grompiex iS the calculated Gibbs free energy of the inhibi- 
tor-fragment-enzyme complex, Greceptor Is the Gibbs free energy for 
the fragment-enzyme complex, and Gigang is the Gibbs free energy 
calculated for the inhibitor compounds. The free energy values for 
each of these terms were estimated as the sum of the four terms: 


G= Eu T Gs sly a Giese — | Snmnede (7) 


in which Ey, is the molecular mechanics energy of the molecule 
expressed as the sum of the internal energy of the molecule plus 
the electrostatics and van der Waals interactions; G,.,,, is the polar 
contribution to the solvation energy of the molecule; G,,,,oi, is the 
nonpolar solvation energy; T is the absolute temperature; and S is 
the vibrational, rotational, and translational entropy of the mole- 
cule estimated by normal mode analysis, using five regularly 
spaced snapshots from the first 100 frames in the last 2 ns of the 
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MD simulations. lonic strength in normal mode calculations was 
set as 0.1M. 


Fragment merging and docking: A fragment merging strategy 
was employed based on the binding conformations observed from 
MM/PBSA studies (Figure 6). The proposed compounds were de- 
signed by merging the structures of lead inhibitors with fragment 
F8. The nitroaromatic ring was deleted, and the remaining pipera- 
zine ring was attached to several different sites on the lead inhibi- 
tors using a variety of different linking groups. With the goal of re- 
taining the lead binding position of the inhibitor scaffold after 
docking the merged hybrids, we investigated several alternative 
linking strategies, including linking to the naphthalene ring, the 
chiral methyl group, and the p-toluidine (scaffold A) or methylene- 
dioxybenzene (scaffold B) by alkyl groups, ester groups, and amide 
groups. The binding conformations of the virtually generated 
hybrid collection were predicted using Glide docking with SP scor- 
ing. The structures of PLpro bound with GRL-0617S (PDB ID: 3E9S) 
and GRL-0667S (PDB ID: 3MJ5) were used in these docking studies 
and prepared using the Schrédinger protein preparation software 
as described above. The docking site was defined as a 10A box 
centered on the original binding site of the lead inhibitors. The 
hybrid compound library was prepared as described above for 
fragment and inhibitor docking. Proposed compounds discussed 
above were selected by score and visual comparison with the crys- 
tallographic position of the lead inhibitor and the predicted bind- 
ing conformation of the fragment. 


Supporting Information 


The replicate plot, ligand efficiency indices, binding free energy 
calculations, and NMR spectra of fragment hit compounds are in- 
cluded in the Supporting Information. 
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